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'sj" I Abstract 

We consider a scalar field (called 0) which is very weakly coupled to thermal bath, 
and study the evolution of its number density. We use the Boltzmann equation derived 
from the Kadanoff-Baym equations, assuming that the degrees of freedom in the ther- 
mal bath are well described as "quasi-particles." When the widths of quasi-particles 
^ . are negligible, the evolution of the number density of (j) is well governed by a simple 

^ I Boltzmann equation, which contains production rates and distribution functions both 

evaluated with dispersion relations of quasi-particles with thermal masses. We pay par- 
ticular attention to the case that dark matter is non-thermally produced by the decay of 
particles in thermal bath, to which the above mentioned formalism is applicable. When 
the effects of thermal bath are properly included, the relic abundance of dark matter 
may change by 0(10 — 100 %) compared to the result without taking account of thermal 
effects. 



1 Introduction 



In particle cosmology, it is inevitable to consider the behavior of quantum fields (or, in other 
words, particles) in thermal bath because the universe was filled with hot plasma in the early 
epoch. The detailed thermal effects depend on how the particle of our interest, which we 
call 0, interacts with degrees of freedom in thermal bath. Importantly, even if is so weakly 
interacting that it is not in thermal equilibrium, there can be non-negligible thermal effects on 
its dynamics. This is because the interaction rate of surrounded by thermal bath depends 
on the properties (in particular, the dispersion relation) of the degrees of freedom in thermal 
bath, which can of course be significantly affected by thermal effects. 

One important example of such a very weakly interacting particle is non-thermal dark 
matter which is produced by the decay of particles in thermal bath. Although the existence 
of dark matter is strongly suggested by various cosmological observations [1], particle-physics 
properties of dark matter, as well as its production mechanism in the early universe, have not 
been understood yet. Various particle physics models including dark matter candidate have 
been proposed so far, like supersymmetric models, universal extra dimension models, and so 
on [2]. In the future, it is hoped that those models are tested by high energy experiments as 
well as by cosmological observations. In particular, the candidate of the dark matter particle 
may be discovered and studied by the LHC experiment as well as future linear colliders, 
based on which a large class of dark matter models are discriminated. For this program, 
precise theoretical calculation of the relic abundance of the dark matter candidate should 
be performed by using information about newly discovered particles [3j. The present dark 
matter density is very accurately determined by the WMAP collaboration as [3]: 

fic/i^ = 0.1126 ± 0.0036, (1.1) 

with h being the Hubble constant in units of 100 km/sec/Mpc, so the dark matter abun- 
dance is now known with 0(1 %) accuracy. Thus, it is desirable to establish methods of 
calculating the dark matter density at the same level of accuracy. For this purpose, detailed 
understanding of thermal effects on the dark matter production process is required. 

In the present study, we pay particular attention to the case that dark matter particle is 
non-thermally produced by the decay of heavier particles in thermal bath. There are many 
examples of such non-thermally produced dark matter, like gravitino [5j, axino [B], a singlet 
field , the right handed sneutrino [S] , and more generally, the recently proposed "freeze-in" 
particles [9]. In such a scenario, the relic density of dark matter is determined at the cosmic 
temperature comparable to the mass of decaying particle and is insensitive to the thermal 
history before that. In the previous studies, the production rate of dark matter has been 
calculated by using the decay rates of particles estimated in vacuum. However, in the actual 
situation, the particles decay in the thermal bath, so the production rate taking account of 
the thermal effects should be properly used in the calculation of the dark matter density. 
As we will see, the thermal effects may significantly change the resultant abundance of dark 
matter. 
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In this paper, we raise the question how important the thermal effects are in the pro- 
duction process of non-thermal dark matter. To answer this question, we first study the 
properties of Boltzmann equation derived from the Kadanoff-Baym equations [10] under the 
assumption that the production of dark matter does not affect the thermal bath. In partic- 
ular, unless the dark matter production is almost kinematically blocked by thermal masses 
during the time when the production of dark matter is most effective, the full Boltzmann 
equation can be reduced to a simplified form which has the same structure as the conventional 
Boltzmann equatioi^J but constructed with "thermal masses" of particles in thermal bath. 
(As we will see, such a simplified Boltzmann equation is obtained by taking the "zero-width 
approximation" of particles in thermal bath.) We evaluate the relic density of dark matter by 
solving (i) the full, (ii) zero-width approximated, and (iii) conventional Boltzmann equations. 
Comparing the three results, we discuss how important the thermal effects are and when the 
zero-width approximation breaks down. We will see that the dark matter abundance may be 
reduced by 0(10 — 100 %) compared with the result of calculation where the thermal effects 
are neglected. 

The organization of this paper is as follows. In Section [21 the relevant formulae to study the 
evolution of the number density of non-equilibrium particles are summarized. In particular, 
properties of the Boltzmann equation to be solved are discussed. Then, in Section |3l we apply 
the formalism to the non-thermal dark matter production process. We numerically solve the 
Boltzmann equation and discuss how important the thermal effects are. Section [His devoted 
to conclusions and discussion. 



2 Formalism 

First, let us introduce the formulae and equations relevant for our analysis. Although many of 
them can be found in literature (see, for instance, [T0l[T2l[l3l[Ta[T5l[l6l[I71[T8l[l9l|20l^ 
we summarize the relevant equations to make this paper self-contained for the sake of readers. 
We assume that the thermal bath have a common temperature T, which have a large degrees 
of freedom, and the back reaction to the thermal bath from the production of is negligible. 

In this paper, we study the evolution of the number density of a scalar field coupled to 
scalar fields in thermal bath, which are denoted as Xi- We introduce the interaction of the 
following form: 

n 

^int = 9(i)Ylxi = 9<pO[xo,Xi,---], (2.1) 
1=1 

where g is a. coupling constant and, for the convenience of the following discussion, we in- 
troduced the operator O = YliXi- The interaction of is assumed to be extremely small. 

In this paper, the "conventional" Boltzmann equation refers to the Boltzmann equation evaluated with 
zero-temperature dispersion relations pT| . 
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i.e., gm"^^^ <^ 1, where is the mass of 0Ej In the following, we study the effects which 
are leading order in g. For simplicity, we consider the case that O is given by a product of 
scalar fields. However, the extension of the formalism to the case that O includes derivatives 
of scalar fields is straightforward. 

The evolution of the number density of (p in the early universe is governed by two effects: 
one is the production of (p due to the decay and scattering processes and the other is the 
cosmic expansion. We discuss these effects separately. 

Because of the weakness of the interaction, (p can be regarded as (almost) free particle!^ 
and the number density operator is given by 

Mt) = 7^ ■■ \kt; ^)kt; -k) + uj'Jit; k)0(t; -k)l :, (2.2) 

where the "dot" denotes the derivative with respect to time, : ■ ■ ■ : is the normal ordering 
and 



c^k = -^k^ + ml. (2.3) 

In addition, 

(2.4, 

where is the field operator for (p and is the volume of the system (where we have adopted 
box normalization). For the later convenience, the number density for each momentum 
eigenstate is defined. Then, the expectation value of the number density of (p is obtained by 
using the density matrix p: 



where, for an operator A, 



iVk(t) = (iVk(t)), (2.5) 



(i) = tr[pi]. (2.6) 



The X sector is in the thermal bath while (p is always out of thermal equilibrium. In 
particular, we are interested in the case that (p is initially absent in the system. Thus, we 
assume that the initial density matrix (at t = ti) is given by the direct product of density 
matrices of two sectors: 



P = P^,i ® Px- (2-7) 

In the case of freeze-in dark matter with a weak scale mass, typically a coupling of O(10~^'^) is necessary 
to obtain the correct dark matter abundance [HI |9] . 

'^•^We assume that the self-interaction of is absent or weak enough to be neglected. We also assume that 
the thermal effects on the expectation value of 4> is negligible, which is the case in the example in Section |3l 
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We set ti = without loss of generality. The x sector is in the thermal bath (with the 
temperature T), so is given by 



Px = e 



(2i 



with being the Hamiltonian for the x sector. On the contrary, j determines the initial 
distribution of (p. We assume that it has translational invariance, i.e., [P,p0^j] = 0, where P 

is the momentum operator. Furthermore, we assume = 0)) = {(f){t = 0)) = 0. 

In order to calculate the evolution of N]i{t), we define the Hadamard propagator and the 
Jordan propagator: 



Gi{t,t';k) ^ (<^(t;k)<^(t';-k)) + (<^(t';-k)<^(t;k)), 
Gtit, t'; k) ^ k)0(t'; -k)) - -k)0(i; k)). 



(2.9) 
(2.10) 



As can be seen from Eqs. (12. 2p and (12. 5p . the expectation value of the number density is given 
by 



iVk(t) 



{dA'+ool)GUt,t';k) 



Cu 



(2.11) 



where Ck is normal-ordering constant. In the weak coupling limit, Ck = |. The Hadamard 
propagator and the Jordan propagator satisfy the following equations, namely the Kadanoff- 
Baym equations [TOl 



{d^ + ul)G'jit,t'-k) 
{d^ + ul) Gi{t,t';k) 



- fdTUt{t-T;k)G^j{T,t';k), 

- fdTUt{t-T;k)Gi{T,t';k) 
Jo 

rt' 

-% / (irn^(t-r;k)Gj(r,t';k), 



(2.12) 



where 

nt(t;k) = -z0(t)(nt(t;k)-nt(t;k) 

<(t;k) = nt(t;k) + nt(t;k). 
Here, at the leading order (i.e, ©(gf^)). 



Ht(t;k) 



nt(t;k) 



tr[e-^x/r] 



-tr 



tr[e-^x/r] 



-tr 



e-^^i^ d{t-k)d{^--k) 
e-^x/r^^O; -k)d(t; k) 



(2.13) 

(2.14) 
(2.15) 

(2.16) 
(2.17) 



with 



dit;k) = L'^/^ j £xe''^-''d{t,x). (2.18) 



Let us define the Fourier transformations 

n^(w,k) = y"cite^'"*n^(t;k), (2.19) 

with Ux = n^et, n^, n<, and n>. Then, from Eq. fl27[D . 

y 27r w — + «0 

By using the relation {u + iO)^^ = Voo^^ — ind^u) (with V denoting the principal value), we 
obtain 

QUUu,k) = -\ [nt(u;,k) -nt(a;,k)] , (2.21) 

and 

nj(u;,k) = -2coth(|,) 9nt(c^,k), (2.22) 

where we have used the so-called Kubo-Martin-Schwinger (KMS) relation n>(ci;,k) = 
exp(a;/T)n<(a;,k) in deriving Eq. (12:221) . 

U sing the fact that the Jordan propagator is time translational invariant within our 
setup^J G'j(t, t'] k) = Gj(t — t', 0; k), it can be expressed in terms of the spectral density 



.#5 



Gj(t; k) = G% 0; k) = 1 ^e--V<^(a;, k), (2.23) 

and the solution to Eq. f l2.12p is given by 

. -2^>ntt(w,k) 

[u;2_^2_3^nt,(u;,k)]2 + [c>nt,(u;,k)P 

We note that the initial condition is given by Gj (0; k) = and S^Gj (0; k) = —i because of 
equal time commutation relations. 



"^■^This is because, under no self interactions of (j) and truncating the perturbative expansion at ©(.g^), the 
spectrum of is determined only by the thermal bath regardless of the number density of (f>. The time 
translational invariance of the Jordan propagator is not a general property. See also [12] . 

"i^^The spectral density should not be confused with the density matrix 



5 



With the initial condition of the Jordan propagator, the Hadamard propagator satisfying 
the Kadanoff-Baym equations fl2.12p and (12.131) is obtained as follows |I2]: 

Gi{t, t'; k) = G'^°'^{t, t'- k)+ [ dh [ dt2G'i{t - h- k)nt(ti - ta; fc)Gj(t2 - t'- k), (2.25) 

JQ Jo 

where the homogeneous solution is given by 

G|°-(t,t';k) = -GUs,s';k) dA'G'^{t;k)G%t';k) 



s,s'=0 



dsGiis, s'; k) {dt + de) Gj(t; k)Gj(t'; k) 



s,s'=0 



-dA,Gt,{s,s';k) G'J(t;k)Gj(t';k). 

s,s'=0 



(2.26) 



In the calculation of the production rate of (p in thermal bath, the most important effect 
is the shift of the pole of the spectral density because its imaginary part gives the production 
rate. Here, we are interested in the case that the interaction of is so small that ^ |n 
Then, the spectral density can be well approximated by the Breit-Wigner form: 



ret I 



(BW) 



'u,k) 



2u;r^(k) 



{co^-coiy + {ur^{k)y 



(2.27) 



where 



r^fk) 



Sntt(a;k,k) 



(2.28) 



Although has the argument k in our expression, it depends only on |k| because of the 
rotational invariance of the thermal bath. Here, we neglected the correction to the real part 
of the pole, which is expected to be irrelevant. In addition, notice that r<^(k) is of 0{g^), and 
is much smaller than Wk- Then, from Eq. (I2.23p . the Jordan propagator (for t > 0) is well 
approximated as 



tG'i{t;k) 



sincjkt r^(k)t/2 



t>o 



Wk 



(2.29) 



resulting in the following expression for the expectation value of the number density defined 
in Eq. (12. 5p . at leading order in F^/cJk, 



(BW) 



it) = /B(u;k)(l-e-^^«0 
1 



Gi{s,s'-k] 



s,s'=0 



ujl+dsds>GUs,s';k] 



2uv 



s,s'=0 



-r^{k)t 



(2.30) 
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with 



fBico) 



(2.31) 



jY^Coii) 



Equivalently, irrespective of the initial condition, one finds 

r^(k)[/B(u;k)-iVk], (2.32) 

from which we can see that T^(k) can be regarded as the production rate of due to the 
decays and scatterings of particles in thermal bath. (Here, the superscript "(Coll)" implies 
that this is the collision term in the Boltzmann equation.) 

As we have mentioned, there is another effect on the evolution of the 0's number density, 
which is the expansion of the universe. Effect of the cosmic expansion can be easily evaluated 
by taking into account the red-shift of the momentum and we obtain 



^(Exp) 



Hk 



dk 



(2.33) 



where the superscript "(Exp)" is for cosmic expansion, and H denotes Hubble parameter of 
the expanding universe. We assume that the energy density of the thermal bath is much 
larger than that of 0, and that the Hubble parameter depends only on the temperature T. 
Combining two effects, the Boltzmann equation to be solved is 

9iVk 



Nu-Hk- 



or, for the total number density 



dk 



T4k;T) [/B(wk;T)-iVk] 



(27r)= 



(2.34) 



(2.35) 



the Boltzmann equation is given by 



dt 



n 



(Coll) 



d^k 
(27r) = 



r^(k;r) [/B(a;k;r)-iVk] 



(2.36) 



In the above equations, we explicitly show that and /b depend on T, which should be iden- 
tified with the cosmic temperature. Thus, the most important quantity to study the evolution 
of the abundance of is the production rate T^{k;T) given in Eq. f l2.28p or, equivalently, 
53Hfgt given in Eq. fl2.2ip . With the operator O given in Eq. (12. ip . the leading contribution 
to '^Hfct is given by 



3Htt(wk,k) = — 



n 



d^Pi 



9_ 

' 2 



n 



i2nr 

d^P: 



G^<ip^,P^ 



Gl'{p^,Pi 



27r)^5(a;k-E.rf)5^'^(k-E.P.) 
(27r)^5(a;k - E.J'°)5^'Hk - E.P. 



(2.37) 
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Here, all the Xj-fields are taken to be independent. If there exist identical fields in {xi} in 
Eq. (12. ip . a symmetry factor is needed. In the above expression, the functions G> and 
are defined as 

G^»(/,p) = ^^^J^^/^^ I c^^xe-P°-°-'P-tr[e-^x/^x.(3:°,x)x.(0,0)], (2.38) 
G^-(pO,k) = tr[e-k/T] / ^'3:e-P°"""'^"tr[e-^>^/^x.(0,0)x.(a:°,x)]. (2.39) 

These functions satisfy the KMS relation G>'(p'^,p) = e^" {p^ , p) , whereas their differ- 
ence, the Jordan propagator, is expressed in terms of the spectral density: 

Gfip', p) = Gl'ip', p) - G^'(/, p) = p,^(p°, p). (2.40) 

Thus, they are given by 

Gl'ip',p) = (/b(/) + 1)Px.(/>p), (2.41) 

Gl'{p',p) = /b(/)p,.(/,p). (2.42) 

If the quasi-particle picture is applicable to Xi, which we assume in the following analysis, 
the spectral density can be approximated by the Breit-Wigner form: 

^(Bw)/ ^ ^ 2p°r^^(p; r) 

(/^-fi2^(p;T))2 + (p0r^^(p;T))2' 

where Q^. and T^. are real and imaginary parts of the pole of the propagator, and are related 
to the real and imaginary parts of the self energy of Xi- Contrary to the case of 0, we need 
to take account of the shirt of the pole for particles which are thermalized. Thus, Q^. may 
significantly deviate from the frequency satisfying the on-shell condition in the vacuum. In a 
large class of models, including the case discussed in the following section, (p;T) can be 
well approximated as 

nl{p;T) = p' + ml{T), (2.44) 

where . (T) is given by the sum of bare and thermal masses. 

The effect of non-vanishing F^. will be numerically studied in the next section. Here, 
we comment that, if F^- is small enough, the evolution of the total number density of is 
governed by a simple differential equation. If the interaction is perturbative, it is usually 
the case that F^.(p;T) ^ m^-(T). In such a case, the (Breit-Wigner) spectral density is 
well approximated as p^^ Xp°,P;T) ~ 27rsign(pO)5(p°^ - - ^x.(^))' ^^^^^ Gl'{u,k) 
and G<'(u;,k) are also (approximately) proportional to the 5-function. We call this limit as 
zero- width limit. Then, the collision term in Eq. fl2.36p becomes 



^(c°ii) 

11' A. 



+ /B(p°))sign(p°) - n /B(P°)sign(p°) 



[/B(fc°)-iVk], (2.45) 
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where 



rfn,.(p.) = |^5(p- -^L(^))' (2-46) 

and dll^{k) is defined in the same way. (Here, we have introduced the four- component vector 
as Pi = {Pi,Pi)-) Notice that, in Eq. f l2.45p . the p° integration is performed in the region 
— oo < < oo, while integration is for > 0. 

It is notable that the collision term given in Eq. fl2.45p includes all the relevant scattering 
and decay processes (and their inverse processes). Regarding Xi scalar particles with 
masses m^^{T), the integrand of the collision term becomes non- vanishing if (and only if) 
the momentum configurations are kinematically allowed. In addition, Eq. (12.451) contains 
the effect of induced emission. For example, if the scattering process Xiipi) ' ' ' Xiipi) ^ 
4'{k)x'i{(li) ■ ■ 'Xpi^p) is kinematically allowed, the collision term contains 



^0 ^ ^7[n^nj^"M[n<^"(^/)]<^°^(^) 

{lujBipmUfi^ + /B(g?))](i + iVk) - [n.(i + /B(p°))][n//B(g?)]iVk}, 

(2.47) 

where we have used the relation /b(— w) = —(1 + /b(u;)). The right-hand side of the above 
equation has the same structure as the collision term in conventional Boltzmann equation; 
however, notice that the thermally corrected dispersion relations should be used in evaluating 
both the phase-spaces and distribution functions. 



3 Application to non-thermal dark matter production 

Now, let us apply the formalism to the non-thermal dark matter production scenario, which 
is recently called freeze-in scenario, regarding as dark matter. In such a scenario, the 
dark matter particle is always out of thermal equilibrium because of the weakness of its 
interaction, and is produced by the decay of particles in thermal bath, Xi- Thus, this is the 
situation where we can safely use the formalism discussed in the previous section. We pay 
particular attention to the question how large the thermal effect can be. For this purpose, we 
numerically calculate the production rate and the relic abundance of (p using the formalism 
presented in the previous section. 

Here, we consider the simplest form of the interaction term, which is 

£int = fl'0C»[xo, Xi] = g<PXQXi- (3.1) 

For simplicity, we concentrate on the case that 0, xoi ^^nd xi ^i-re all real scalars, with masses 
satisfying m^^ > + m^^. In our numerical calculations, we take m^^ = 100 GeV and 
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mj(,j = 0. Then, in the absence of thermal effects, can be produced by the two-body decay 
process Xo ~^ As we will see in the following, such a decay process may be blocked, or 
even the "decay" process xi — ^ (pXo may occur, once thermal effects are taken into account. 

At the leading-order in g, is produced only by 1 <H- 2 "decay" processes. Using the 
formulae given in the previous section, we obtain 

9^ f d'^Q 

r^{k;T) = —J ^^[l + /B(go) + /B(wk-go)]Pxo(?o,q)Pxi(^k-go,k-q). (3.2) 

As we have discussed in the previous section, when the quasi-particle picture is applicable, 
the thermal effects on Xi are imprinted in Q^- and T^.. Importantly, these quantities depend 
on how Xi interacts in thermal bath. To make our discussion concrete, we adopt the following 
form of the interaction 
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Ant = -^ai^xt- (3.3) 



Then, at the leading order in g^^, the thermal mass of Xi is obtained as 

= mz + / - — 3.4 

^' 2 J (27r)3 ' ^ ^ 



where u^' = ^/m^^ -|- p^. In the following analysis, we consider relatively large value of g^-^. 
Thus, we include the NLO contribution to the thermal mass of Xi- In the present analysis, 
as we mentioned, we take m^-^ = 0. Then, the thermal mass of Xi is given by, at NLO in g^-^ 

m 



3 9x1 



(3.5) 

/24; ^ ' 

Note that we drop the momentum dependence, since, at this order, the only contribution is 
a tad pole diagram. In addition, F^. is given in the integral form as [15j 



TT 



4 



^^^P' ^ 2fi,,(p)^^»y (27r)32fi,^,i y (271)32^,^,2 7 (27r)32fi,,,3 

(27r)35(p + - q2 - q3)5(fi^^(p) + fi^^^^ - fi^^ 2 - ^^x^.s) 

[/b,i(1 + /b,2)(1 + /b,3) - (1 + /b,i)/b,2/b,3] , (3.6) 

where, for the simplicity of the equation, we defined ^xi,J = ^Xii^j) /b,J = /B(^Xi.-') 
(with J= 1-3). 

With the above formulae, we can follow the evolution of the number density of in the 
early universe and calculate the relic abundance. For this purpose, it is convenient to define 
the "yield variable" as 

Y, = ^, (3.7) 

s 



Although Xi couples to other particles in the thermal bath, here we assume that the dominant effect on 
the dispersion relation of Xi is from its self interaction. 
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where s = ^g^T^ is the entropy density. (Here, is the effective number of relativistic 
degrees of freedom; in our numerical calculation, we use (7* = 100.) We are interested in the 
case that A^k <^ fsi^^k) and, in such a case, we can neglect the term proportional to A^k in 
the right-hand side of Eq. fl2.36p . Then, the present value of is given by 



d log z 



dlogz 



where 



dlogz 



n 



(Coll) 



sH 



d^k r^(k;T)/BM 
(27r)3 sH 



with 



m 



z = 



XO 



T 



(3.9) 



(3.10) 



Notice that the variable z increases as the universe expands. We will show how dY^/ dlogz 
and Y^ behave as a function of the cosmic temperature T in the following. 

Before showing the numerical results, it is instructive to consider the zero-width limit; in 
the present case, Eq. f l2.45p becomes 



n 



(Coll) 



d^p 



m 



XO 



(27r)3fi^„(p 



■rxo->0Xi/B(^Xo(p))[l + /xi(P)] 



m 



XI 



-fxi^<^Xo/B(^Xi(P))[l + /xo(P)] 



(3.11) 



and r 



Xl-5>(/>X0 



where we have used the relation A^k ^ /B(i^k)- In the above expression, T^^^^^^ 
are decay rates of Xo and Xi calculated with the conventional Feynman rules but with ther- 
mally corrected dispersion relation; for example, T^^^^^-^ is given by 



xo->-<t>xi 



9' 



IGvrm^g 



V ^xo ^xo 



(3.12) 



In addition, /^^ and /^^ are "averaged" distribution functions of Xo and Xi in thermal bath, 
respectively. For example, /^^ is given by 



/xo(p) 



rfcos 6'cm/b(^xi(IpI'^cm)), 



(3.13) 



where E^j^{\p\,9cm) denotes the energy of Xi emitted from Xo (carrying momentum |p|) to 
the direction ^cm relative to p (with 9cm being defined in the rest frame of Xo)- Notice 
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dY^/cllnz dY^/dlnz 




0.001 0.01 01 1. 10. OOOl 0.01 0.1 1. 10. 



Figure 1: dY^p/dlog z as a function of z = m^^g/T for m^/m^^ = 0.95, g^^ = 0.1, g/my,^ = 10~^^, 
and g^ = 100. The left figure is computed with g-^^^ = 0.3 and the right one is with gy,-^ = 2. The 
red sohd hues are the results of full calculation while the black dashed ones are zero-width results. 
Results of conventional Boltzmann equation are also shown in blue dot-dashed lines. 



that, in the calculation of Ey-^{\p\,9cM), dispersion relations including the thermal masses 
should be used. Eq. (13. lip indicates that the production process of (p is active only when 



fhy^{T) > + my^{T) or m^j(T) > + rhy^^{T)]^ otherwise, the production of is 
kinematically suppressed due to the thermal mass. As we will see below, this is indeed the 
case. 

Now we are at the position to numerically calculate the production rate of in thermal 
bath. First, let us discuss how dY(j,/d\ogz given in Eq. (13. 9p behaves. In Fig. [H we plot 
dY^/d\ogz as a function oi z = m^^^jT . Here we take m^/rriy^^ = 0.95, g^^^ = 0.1, and 
Qy-^ = 0.3 (left) and g^-^ = 2 (right). In each figure, the result of the full calculation, which 
is calculated by using Eq. (13. 9p including the effects of the widths of Xi, Txi' shown in the 
red solid line. In the same figure, in the black dashed line, we show the result taking the 
zero-width approximation (see Eq. (I3.1ip ). The result of the calculation with conventional 
Boltzmann equation (which corresponds to the zero- width approximation with rh-^^ — )• rriy.) 
is also shown in the blue dot-dashed line. 

As one can see, dV^/dlogz is significantly suppressed when 0.027 < m^^/T < 1.2 (for 
g^-^ = 0.3) or 0.16 < rriyg/T < 6.3 (for g^-^ = 2). This is due to the fact that, at leading 
order in F^./m^., the production process of is kinematically blocked because of the thermal 
mass; the upper and lower edges correspond to the temperature where m^^ ~ + rh^g and 
~ + rhy^ are realized, respectively. In the suppressed region between the edges, the 
production is from the off-shell effects of Xi (especially Xi); namely due to non- vanishing 

mm- 

It is notable that, in the region where there is no kinematical suppression, the result of 
zero-width approximation agrees remarkably well with that of full calculation. Thus, the 



'^''In the present case, (j> should play the role of dark matter, and hence is stable. Thus, we do not have to 
consider the case that > m^g (T) + m^^ (T) > m^g + m^g . 
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effect of T^. can be safely neglected in the calculation of the relic abundance of 0, unless 
the production channel is almost kinematically blocked at the epoch when the production 
of (j) is most active. If the zero-width approximation can be adopted, the calculation of the 
production rate is extremely simplified. In the opposite case where (j) is mostly produced 
at the epoch of kinematical suppression, production processes due to T^. ^ may not be 
negligible in the calculation of the relic abundance of (p. A typical example is given in the 
right hand side of Fig. [H 

From Fig. [H we can also see that, if we perform a naive calculation without properly 
taking account of the thermal effect, the production rate may significantly deviate from the 
correct value. In particular, even at > + rh^g or rh^^ > + m^^ , the discrepancy 
between the results of full and conventional calculations is sizable in some parameter region. 
This is because, at m^^ > + m^^, the dominant production process of is Xi 0Xo 
contrary to the case of low enough temperature, and at > + m^i' phase space 
becomes smaller than that in the case of conventional calculation due to the thermal mass of 
Xi- 

In the non-thermal dark matter production due to the decay, the dominant production 
occurs when m^^jT 1—5; for such a temperature, the actual value of dY^/d log z may receive 
sizable thermal effects. Thus, for the accurate calculation of the relic density of non-thermally 
produced dark matter, it is dangerous to neglect the thermal effects. In Figs. [2] and [HI we 
show the evolution of the yield variable Y^. In the same figure, we also show the result of the 
naive calculation neglecting the thermal effects. As one can see, the evolution of Y^ changes 
once the thermal effects are taken into account, and the resultant value of Y^ is suppressed 
compared to the result obtained by neglecting the thermal effects. In particular, for the case 
of (yf^j = 2, which gives relatively large value of F^^, significant amount of is produced even 
during the period of kinematical suppression and the resultant Y^ differs from the zero-width 
result. This is because the production channel is suppressed during the time m^^/T ~ 1 — 5, 
where the most effective production could occur if there were no thermal effects. Therefore, 
the production due to the off-shell effects with F^. ^ becomes important. In such a case, 
the resultant value of Y^ does not agree with the result of zero-width approximation nor that 
of the conventional Boltzmann equation. 

In Fig. m we plot the ratio of the yield variable calculated with and without thermal 
effects, 



R ^" 



(conventional Boltzmann equation) 



(3.14) 



as a function of g^-^^ and nifp/m^g, using the zero-width approximation. In the present set up, 
the present number density of decreases by taking account of the thermal effects. As one 
can see, Y^ is more suppressed when g^^ becomes larger or when the mass of becomes closer 
to that of xo- 

Before closing this section, we note here that there is another possible mechanism of 
enhancing the production rate. With the emission of a massless particle in the thermal bath 
from the initial or final state particles (corresponding to Xo or Xi in the present setup), the 
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5 10 15 20 0.01 0.1 1. 10. 

Figure 2: The evolution of the yield variable for m^/m-x^^ = 0.95, g-^^ = 0.1, g-^^ = 0.3, 
d/'^xo ~ 10~^^) aiid gi* = 100. The red solid line is computed with T^^ ^ and the black dashed 
one is done with F^. = 0. The blue dot-dashed one represents the conventional result (without 
thermal effects). 




5 10 15 20 0.1 1. 10. 



Figure 3: Same as Fig. [2] except for g-^^ = 2. Notice that, in the left figure, the line for the 
conventional result is almost parallel to the vertical axis. 



production rate of the non-thermal dark matter at 2; < 1 may be affected, and may become 
larger by the factor of a few [23]. In realistic models of non-thermally produced dark matter, 
the particles corresponding to xo and xi niay couple to massless gauge bosons, which may 
enhance the production rate at 2; < 1. In the present model, however, there is no massless 
particle responsible for such an enhancement. In addition, for the calculation of the relic 
density of the non-thermally produced dark matter, such an effect may not be important 
because, as we have shown, the relic abundance of the non-thermally produced dark matter 
is determined at z ~ 0(1). 



14 



R 




0.2 0.4 0.6 0.8 1.0 



Figure 4: The contour plot of i? as a function of g^^^ and m^/my^^, taking =0.1. 

4 Conclusions and Discussion 

In this paper, we have discussed the evolution of the number density of a particle which 
is coupled to thermal bath very weakly and hence is in non-equilibrium. We first solved 
the Kadanoff-Baym equations for the case that (i) the effects of (j) on the thermal bath is 
negligible and (ii) the self interaction of (f) is sufficiently weak. Then, we derived Boltzmann 
equation describing the evolution of the 0's number density for the case that (iii) the real 
part of the 0's self energy is (almost) unchanged by thermal effects. We then studied the 
properties of Boltzmann equation, assuming that the degrees of freedom in the thermal bath 
are well described as the quasi-particles. In particular, in the situation that the widths of 
quasi-particlcs arc negligible, the evolution of the number density is well described by the 
Boltzmann equation in the familiar form, which contains the matrix elements, phase-space 
integrals, and distribution functions evaluated with the "on-shell" condition modified by the 
thermal effects. 

Then, we have applied the formalism to the scenario in which dark matter is non-thermally 
produced from the decay of "particles" in thermal bath, regarding (j) as the dark matter. 
In such a case, the above conditions (i), (ii), and (iii) are satisfied, and the formalism we 
have studied can be safely apphed. We calculated the number density of the non-thermally 
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produced dark matter, taking account of the effects of thermal bath. Because of the change 
of the dispersion relations of "particles" in thermal bath, the production rate of the dark 
matter may significantly change. In particular, in some cases, the decay process to produce 
0, which is kinematically allowed in the vacuum, may be blocked because the particles in 
thermal bath acquire thermal masses, which changes the mass relation among parent and 
daughter quasi-particles. Numerically, we found that, if the proper Boltzmann equation with 
the thermal effects is used, the dark matter density may change by 0(10 — 100 %) compared 
to the results of calculations neglecting the change of the dispersion relation of the "particles" 
in thermal bath. We have also studied the effect of the width in the spectral density of the 
quasi-particles in thermal bath, and found that the zero-width limit can be safely taken if 
the thermal blocking of the production of does not occur at the time when the production 
of is most effective. 
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